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Abstract. We report the firm discovery of solar-like oscillations in a giant star. We monitored the star £ Hya (G7I1I) 
continuously during one month with the CORALIE spectrograph attached to the 1.2m Swiss Euler telescope. The 
433 high-precision radial-velocity measurements clearly reveal multiple oscillation frequencies in the range 50 - 
130 fiHz, corresponding to periods between 2.0 and 5.5 hours. The amplitudes of the strongest modes are slightly 
smaller than 2 ms" . Current model calculations are compatible with the detected modes. 
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1. Introduction 

Doppler studies with high-precision instruments and re- 
duction algorithms have been refined dramatically, mainly 
in the framework of the search for exoplanets. These re- 
finements have led to a breakthrough in the observations 
of solar- type oscillations, which have n ow been found re- 
peatedly (Procyon, Martic et al. 1999j (3 H yi, Bedding et 
al. |200l|; a CenA, Bouchy & Carrier |200i]; 5 Eri, Carrier 



et al. |2002|) . The signal-to- noise ratio (S/N) in the oscilla- 



tion frequency spectra is, for each of these cases, so good 
that the resemblance with the solar oscillation spectrum 
is obvious. 

Observations of solar-like oscillations in the giant star 
aUMa have been claimed by Buzasi et al. ( [2000 ), based 
upon space photometry gathered with the star tracker of 
the WIRE satellite. However, the interpretation of these 
reported oscillations frequencies is not straightforward. 
Guenther et al. (2000) find a possible solution in terms 
of a sequence of radial modes with a few missing orders 
for a star of 4.0-4.5 M Q . The interpretation is not sup- 
ported by theoretical calculations by Dziembowski et al. 



(2001). Velocity observations of Arcturus provide evidence 
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for solar-type oscillations with periods from 1.7 to 8.3 
days and a frequency separation of evenly spaced modes of 
Av ~ 1.2 fiRz (Merline, 1999). WIRE data (Retter et al. 
[2002 ), however, points to an excess power at 4.1 /iHz and 
a frequency spacing of Av = 0.8 /iHz. 

In this Letter, we provide clear evidence for the pres- 
ence of solar-type oscillations in the giant star £ Hya 
(my = 3.54). This star has a mass close to M — 3.0M Q , 
and is thus considerably heavier than the Sun. Moreover, 
its luminosity amounts to L ~ 61 L Q and its effective tem- 
perature T G ff — 5000 K, which places the star among the 
giants. In the current Letter we present the first results 
of our study. Detailed modelling will be presented, when 
completed, in a subsequent paper. 

2. The radial-velocity measurements 

The £ Hya observations were made during one full month 
(2002 February 18 - March 18) with Coralie, the high- 
resolution fiber-fed echelle spectrograph (Queloz et al. 
[2001 ) mounted on the 1.2-m Swiss telescope at La Silla 
(ESO, Chile). During the stellar exposures, the spectrum 
of a thorium lamp carried by a second fiber is simultane- 
ously recorded in order to monitor the spectrograph's sta- 
bility and thus to obtain high-precision velocity measure- 
ments. The spectra were extracted at the telescope, using 
the INTER-TACOS (INTERpreter for the Treatment, the 
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Fig. 1. Radial- velocity measurements of £ Hya. The dis- 
persion reaches 4.6ms" 1 . a) The whole time string; b) 
zoom of a small part of the time series. 
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Fig. 2. Power spectrum for £ Hya. Modes are clearly 
present in the range 50-130 /iHz, corresponding to pe- 
riods of 2.0-5.5 hours. Panel a) is the power spectrum of 
the time series using one single spectrum as the reference. 
Panel b) is for the time series derived using one reference 
spectrum per night. The insert in Panel b) shows the win- 
dow function. The frequency scale is made similar to the 
main plot to facilitate comparison. 



Analysis and the Correlation of Spectra) software pack- 
age developed by D. Queloz and L. Weber at the Geneva 
Observatory (Baranne et al. 1996| ). The wavelength cov- 
erage of these spectra is 3875-6820 A, recorded on 68 or- 
ders. By taking about 2 measurements every hour, a total 
of 433 optical spectra was collected. The exposure times 
were typically 180 s and the S/N ratio for all spectra was 
in the range of 110-230 at 550 nm. 

By the use of the optimum- weight procedure (Bouchy 
et al. 2001 ), radial velocities are computed for each night 
relative to the highest S/N spectrum obtained in the mid- 
dle of the night. This method requires a Doppler shift that 
remains small compared to the line-width (smaller than 
100 ms" 1 ) (Bouchy fc Carrier pOOl]) . Since the Earth's mo- 
tion can introduce a Doppler shift larger than 800ms _1 
during a whole night, each spectrum is first corrected for 
the Earth's motion before deriving the radial velocities. 
This is achieved by shifting all spectra with the Earth's 
velocity at the time of observation and by subsequent re- 
binning, so that the spectra all have the same wavelength 
values. From each rebinned spectrum a velocity is derived. 
The mean radial velocity for each night is then subtracted. 
The resulting velocities are shown in Fig. 1. The rms scat- 
ter of the time series is 4.6 ms" 1 and is largely due to 
the oscillations. The mean error on each measurement is 
2.3ms" 1 . 



3. Time-series analysis 

The power spectrum of the 433 data points is presented in 
Fig. |[ In order to show that the excess power in the lower 
panel in the range 50-130 /iHz is not due to elimination 
of power at low frequencies by the reduction procedure, 
we show a power spectrum in the upper panel where no 
correction for drift has been applied. Although more noisy 
the increase in power between 50 /iHz and 130 /iHz is also 
evident in the upper panel. A set of oscillation modes is 
clearly present, with a power distribution that is remark- 
ably similar to that seen in the Sun and other stars on or 
near the main sequence, although obviously at much lower 
frequency. We will show the spectrum displays the char- 
acteristic near-uniform spacing of the dominant peaks. 

To characterize this pattern we have calculated the au- 
tocorrelation of the power spectrum. In order to eliminate 
the effect of the noise, we have ignored all points with am- 
plitudes below a given threshhold (1.2 ms" 1 ). In Fig. | 
the alias at 1 c/d stands out clearly, but in addition a 
spacing Av ~ 6.8 /iHz is present in the power spectrum. 
This is consistent with the visual impression of a regular 
pattern present in the the power spectrum (Fig. |J). 

The expected velocity amplitude for solar-like oscil- 
lations scales as L/M according to Kjeldsen & Bedding 
(1995). Using the stellar parameters from Sect. 4 and a 
mass of 3.0 M Q (Sect. 4) we find w sc,pred = 4.6ms" 1 . The 



S. Frandsen et al.: Detection of Solar-like Oscillations in the G7 Giant Star £Hya 



3 



0.004 r 



0.003 



0.002 



0.001 



o.ooo LL 







5 10 
Freq//xHz 



1 b 



Fig. 3. Autocorrelation of the power spectrum for £ Hya. 
All points with amplitude v < 1.2 ms have been set to 
zero. The dash-dot-dot-dot vertical lines mark the position 
of the 11.574 /iHz daily alias and half of it, and the dash- 
dot line indicates the peak pointing to a frequency spacing 
of 6.8 [iHz. 



observed amplitudes in Fig. g are only about one third of 
this prediction. Recent calculations by Houdek & Gough 
(2002), however, indi cate t hat the simple scaling law of 
Kjeldsen & Bedding ( 1995 ) does indeed not apply. They 
predict a velocity amplitude v — 2.1 ms , slightly higher 
than the ones observed by us for £ Hya. 

The stochastic nature of solar-like oscillations implies 
that a timestring of radial velocities cannot be expected 
to be a set of coherent oscillations and can therefore not 
be reproduced perfectly by a sum of sinusoidal terms. As a 
starting point it is nevertheless a good assumption to try 
to fit the radial- velocity data of £ Hya by such a set of func- 
tions, as the lifetimes are expected to exceed the length 
of the observing run (Houdek & Gough 2002). We have 
performed an iterative fit u sing different methods, among 
which Pcriod98 (Sperl 1998 ) and a procedure by Frandsen 
et al. ( 2001 ). An oscillation with a S/N above 4 is detected 
at 9 frequencies. When a fit based upon these 9 frequen- 
cies is removed from the time series, 5 additional peaks 
still occur but with a too low S /N to accept them without 
additional confirmation (hence we do not list them) . After 
removing the 14 frequencies, only a noise spectrum is left. 
The results for the 9 frequencies are presented in Tabic |l|, 
where the S/N indicated is calculated from the remain- 
ing noise in the amplitude spectrum at the position of 
each mode. The noise is slightly higher in the range of the 
modes than at high frequencies, where a = 0.2 ms -1 . The 
frequencies of the modes with amplitudes above 1.5 mr 1 , 
i.e., of the five highest-amplitude modes, are unambigu- 
ous. Dividing the dataset in two, four out of five modes 
with S/N > 5.0 are present in each set. For lower S/N 
the alias problems lead to a risk that false detections are 
made. Modes with S/N < 5 must be considered with some 
caution. Confirmation of these frequencies is needed by 
additional observations. What is stated above has been 
verified by the analysis of several sets of simulated data 
assuming a variety of lifetimes in order to check the valid- 



Table 1. Oscillation frequencies detected in the radial 
velocities of £ Hya ordered by amplitude. Formal errors of 
the fit are given in parenthesis as errors of the last digits. 
The S/N is calculated by dividing the amplitude by the 
noise in a 10 //Hz bin in the residual spectrum centered on 
the frequency of the mode. Frequencies in italic indicate 
modes where the frequency detection methods disagree 
on the correct alias. The last column gives the difference 
between the observed frequency and the frequency derived 
from Eq. (1). 



ID 


Frequency 


Amplitude 
ms~ 


S/N 


Av 




c/d 


/iHz 




/iHz 


Fl 


5.1344(26) 


59.43 


1.85(23) 


6.6 


0.77 


F2 


6.8366(27) 


79.13 


1.84(23) 


5.8 


-0.86 


F3 


7.4265(29) 


85.96 


1.76(23) 


5.3 


-1.14 


F4 


8.2318(32) 


95.28 


1.65(23) 


5.1 


1.07 


F5 


9.3507(33) 


108.22 


1.59(23) 


6.0 


-0.21 


F6 


8.7399(36) 


101.16 


1.36(22) 


4.5 


-0.16 


F7 


10.0287(43) 


116.07 


1.25(23) 


5.0 


0.53 


F8 


9.0831 (U) 


105.13 


1.24(24) 


4.3 




F9 


8.5339(40) 


98.77 


1.23(23) 


4.1 





ity of the identified modes. The details of such simulations 
will be published in a subsequent paper. 

The first seven modes can be ordered in a sequence of 
modes, which fits the straight line 



v(n) = 1.78 + 7.11(±0.14) 



(1) 



where n is an integer (the order of the mode). Some n 
values are missing. The maximum deviation from the line 
is 1.14 /tHz (Table 1). The regularity seen is similar to the 
results reported for Ar cturus (Merline, 1999) and a UMa 
(Guenther et al. 



2000) 



The first seven modes could all be 
radial modes, although we cannot rule out the possibility 
of alternating degree i — and I = 1 modes. This, how- 
ever, is not consistent with the model presented in Sect. 
4. 

From the present data, we cannot firmly exclude that 
frequency F8 in Table 1 corresponds to the same mode as 
F7. The two modes are resolved, but if the damping time 
is short, they might be different realizations of the same 
mode. 

4. First interpretation 

In order to study the nature of the oscillations detected in 
£ Hya, it is necessary to compare the observed frequency 
spectrum with model predictions, taking into account the 
constraints on the three observational stellar parameters 
T c ff, L and Z . We have redetermined the atmospheric pa- 
rameters and use T off = 5000 ±100 K, L = 61.1 ±6.2 L , 
logff = 2.85, [M/H] = 0.06 ±0.07 (Z = 0.019 ±0.006) and 
vsini = 1.8±1.0 kms -1 . Details of how these values were 
obtained will be reported in a subsequent paper. 

Using the evolution code of Christensen-Dalsgaard 
(1982), evolutionary tracks were produced, spanning the 
error box defined by the uncertainties of T e g, L and Z. 



4 



S. Frandsen et al.: Detection of Solar-like Oscillations in the G7 Giant Star £Hya 




8.0X1CT 
T.„ (K) 



Fig. 4. Hertzsprung-Russell diagram with evolutionary 
tracks for Z = 0.019 for masses 3.07M© (solid line), 2.93 
and 3.15M Q (lower and upper dashed lines, respectively). 
The location of £ Hya in the diagram is given by the error 
box. 



The model tracks were computed using the EFF equation 
of state (Eggleton et al. 1973), OPAL opacities (Iglesias, 
Rogers k Wilson 1992), Bahcall & Pinsonneault (1992) 
nuclear cross sections, and the mixing-length formalism 
(MLT) for convection. 

The evolutionary track passing through the observed 
(T c ffj L) corresponds to a mass of 3.07M Q for Z — 0.019 
(Fig. |]). Oscillation frequencies were calculated for the 
model in that track closest to the observed location of 
£ Hya in the HR diagram. The average separation between 
radial modes in the range 50-100 /TFz is Av t h co = 7.2 /iHz 
in agreement with the observational value. The model fre- 
quencies fit a linear relation for orders 5 < n < 17 or 
50 < v < 125 /kHz given by v = 7.15 + 7.14 n with ab- 
solute values only 1-2 (iHz from the observed frequencies. 
The maximum deviation of the model frequencies from 
the line is ~ 0.5 /.iHz. 

Radial modes are expected to dominate the spectrum 
for giant stars (cf. Dziembowski et al. 2001, Fig. 2). 



Further analysis of the spectrum is beyond the scope of 
the current discovery paper and will be done in a forth- 
coming paper, dealing in detail with the issue of modelling 
and considering also the possibility that £ Hya could be a 
core helium burning star with a smaller mass. 



5. Conclusions 

The main conclusions of this study are as follows. Solar- 
like oscillations have been firmly discovered in the bright 
G7III star £ Hya. The amplitudes of the strongest modes 
are somewhat below 2ms . The observed frequency dis- 
tribution of the modes detected (Table 1) is in agreement 
with theoretically calculated frequencies both in terms 
of the spacing and the absolute values. The modes with 
the largest amplitudes can be well matched with radial 



modes that have an almost equidistant separation around 
7.1 /iHz. 

A most important and exciting result of our study is 
the confirmation of the possibility, suggested by the results 
reported on a UMa and Arcturus, to observe solar-like 
oscillations in stars on the red giant branch. This opens 
the red part of the HR diagram for detailed seismic studies. 
The latter require an accuracy within the range of current 
and future ground-based instruments. Such future studies 
will only be successful if an extremely high stability of 
the instrument is achieved and if one performs multisitc 
observing campaigns covering several months in order to 
resolve the frequency spectrum of the oscillations and to 
eliminate the aliasing problems. 
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